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Abstract 

In a preceding study in the heavy quark limit of QCD, it has been demon- 
strated that the best lower bound on the curvature of the Isgur-Wise function 
£(w) is £"(1) > i[4p 2 +3( / o 2 ) 2 ] > i|. The quadratic term (p 2 ) 2 is dominant 
in a non-relativistic expansion in the light quark, both £"(1) and (p 2 ) 2 scaling 
like (R 2 m 2 ) 2 , where m q is the light quark mass and R the bound state radius. 
The non-relativistic limit is thus a good guide-line in the study of the shape 
of £(w). In the present paper we obtain similar bounds on all the deriva- 
tives of £nr(vj), the IW function with the light quark non-relativistic, and 
we demonstrate that these bounds are optimal. Our general method is based 
on the positivity of matrices of moments of the ground state wave function, 
that allows to bound the n-th derivative £pjji(w) in terms of the m-th ones 
(m < n) . We show that the method can be generalized to the true Isgur Wise 
function of QCD £(w). 
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1 Introduction. 



Using the OPE in the heavy quark limit of QCD, new Bjorken-like sum rules 
(SR) have been obtained |TJ 12 El- It has been shown that the Isgur-Wise (IW) 
function C,(w) is an alternate series in powers of (w — 1), and lower bounds have 
been found on the absolute magnitude of its derivatives. Important ingredients in 
the derivation of the SR are the consideration, following Uraltsev jl] of the non- 
forward amplitude, plus the systematic use of boundary conditions that ensure that 
only a finite number of j p intermediate states (with their tower of radial excitations) 
contribute. 

In particular, it has been found that the n-th derivative is bounded by the (n— 1)- 
th one |2] 

(_i n (n) (1) > ?n+l ( _ ir x^-i) (1) > (1) 

where the second inequality follows from the recursive character of the first one, and 
generalizes the inequality for the slope p 2 = — £'(1) : 

p 2 a \ CO 

that follows from Bjorken [3] and Uraltsev jl] SR. The first inequalities (JTJ) read, for 
the curvature a 2 = £"(1) 

In |3] one has obtained, from a wider class of SR, the following better bound on 
the curvature 

a 2 >^[4p 2 + 3(p 2 ) 2 ] >^ (4) 

where the absolute lower bound (independent of p 2 ) follows from (J2J). Radiative 
corrections to the bounds (JHJ) and (jlj) have been computed by M. Dorsten jB]. 

It has been underlined in that the quadratic term in (@J) has a clear physical 
interpretation, as it is leading in a non-relativistic (NR) expansion in the mass of 
the light quark : 
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> ^(/) 2 (5) 

where £nr{w) denotes the Isgur-Wise function with the light quark in the NR limit. 

It is clear that it is very important to have rigorous bounds on the derivatives of 
the IW function. The main general reason is that the shape of the latter is linked to 
the determination of the CKM matrix element \V c b\ through the exclusive processes 
B -> D^£u. As pointed out in [3], a more quantitative reason is that, beyond the 
first derivative, higher derivatives will play a non-negligible role at the edge of the 
phase space, at high w, in the region where the data are presently rather precise, 
and will become more and more precise in the near future. 

Therefore, it is not only of an academic interest to find bounds on higher deriva- 
tives. To this aim, one can begin by using systematically for higher derivatives the 
method exposed in It is then possible to find bounds of the form (|3J) for higher 
derivatives [Zj. 

However, we have realized that, studying the NR situation, a more powerful 
method can be developed that leads to better bounds, exposed below. With this 
method, one finds for the second derivative, in the NR limit, the bound (jSJ, but 
for higher derivatives £jvi?(l) one finds more complicated bounds involving the m-th 
ones {m < n). It is important to notice that, unlike the method that we have used 
in QCD pQ-jBI, where the IW functions to excited states played a crucial role, in the 
present paper we note that a simple general property (positivity) involving only the 
elastic IW function allows to deduce all the bounds on its derivatives at zero recoil. 

In QCD, in the heavy quark limit one can hope to obtain bounds such that, in 
the NR limit for the light quark, reduce to these better bounds. For the moment, the 
aim of the present paper is to obtain bounds for the derivatives of £nr(w), opening 
the way to study the more complex case of the actual IW function £(w) in QCD. 

In Section 2 we set the relation between the derivatives and the moments 

< 0|,z m |0 >, where |0 > is the ground state radial wave function. In Section 3 
we study the general constraints on the moments. In Section 4 we shift to the 
corresponding constraints on the derivatives and the resulting bounds. In Section 
5 we illustrate our bounds by particularizing to simple potentials and in Section 6 
we conclude. In Appendix A we demonstrate a mathematical identity used in the 



3 



text, in Appendix B we show the frame dependence of subleading moments and 
in Appendix C we demonstrate the optimality of the bounds. In Appendix D we 
deduce a simple formula giving a weaker but completely explicit bound for all the 
even derivatives ^ 2n ^(l), that we compare with our optimal bounds. 

2 Isgur-Wise function in the NR limit : relation 
between its derivatives and moments. 

2.1 Universal NR form factors in heavy- heavy transitions. 

As it is well known from nuclear physics, the first-quantized non-relativistic 
operator corresponding to the electromagnetic current eqQ(x)'~f Q(x) of field the- 
ory is given by the expression eg<5(x — tq), where vq is the position of the active 
quark (see for example [S] ) . Here we are interested in the form factor corresponding 
to a general current Q'{x)YQ{x) for transitions between bound states of the type 
{Qi q) ~ > {Q' i q) with unequal masses. For simplicity, let us begin with the matrix 
element of the fourth component of the vector current J°(0) = Q '(0)j°Q[0). 

Let us write the corresponding non-relativistic (NR) form factor in a general 
frame. The form factor for the transition (Q, q) — > (Q', q) will be given simply by 
the matrix element of the operator 5(tq) : 

F(P',P) =< P'|5(r Q )|P >= J ¥ P ,(r q ,0)*¥ P (r g ,0)dr q (6) 

where yq and r q are respectively the positions of the active and the spectator quarks. 

There is a simple argument to find this non-relativistic expression of the current 
J° = Q'^Q. Its matrix elements between one-particle states of given momenta 
p and p' in the non-relativistic limit, where p'/ m Q' an d p/^q are small, are 
< p',s'|J |p, s >= 5 s r yS . As readily verified, these matrix elements are precisely 
those of the multiplication operator (in configuration space) by the function S(r). 
Indeed, one has : 

<p>,s'\6(r)\p,s> =6 S ,, S J e~^'- r 6(r) e^ r dr = 6 S , S . (7) 
The non-relativistic limit of the current is therefore J° = S(r). 
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The wave functions, that factorize in center-of-mass and internal wave functions, 
are given by 

*p( r «> r Q) = ex P ( m Q r q + m Q r Q) ■ v ] i> 1 ( r ? - r o) 

*P'( r </> r Q') = eX P [* ( m q T q + m Q' r Q') ■ V '] ( r <? - r Q') (8) 



where 



P' 



(9) 

m q + m,Q m q + 

are the non-relativistic velocities and one gets for the form factor (JOJ) the following 
expression, exhibiting Galilean invariance, since the NR form factor F(P, P') is a 
function of the variable (v — v') 2 : 

F(P', P) = |exp [m 5 (v - v') • r,]| tf) = /((v - v') 2 ) . (10) 

Up to now, equations (|B|)-()1L)|) are valid for the current J°(0) = Q'(0)7°Q(0) and 
for any value of the masses rriQ, mcy, m q . If we now assume the hierarchy 

1/R < m q < m Q) m Q i (11) 

we will be in the situation of a heavy-heavy Q — >■ quark transition with a light 
spectator, non-relativistic quark g. The first condition 1/R <C rn. g , where i? is the 
radius of the bound state, ensures that the quark q is non-relativistic, while from 
the second condition m q <C mQ,mQ>, the quark g is light relatively to the active 
quarks Q and Q'. This latter condition implies a heavy quark symmetry SU(2Nh) 
where is the number quarks that are heavy relatively to the quark q. Therefore, 
these conditions imply Isgur-Wise scaling, with all form factors being given by the 
universal form factor /((v — v') 2 ) IjlOj) . Notably, the flavor independence is due to 
the fact that the internal wave function, ipf and ipi become independent of the heavy 
quark mass. 

Expanding the form factor /((v — v') 2 ) in powers of (v — v') 2 we can write then 
/((v-v') 2 )=E i/W(0)(v-v') 2 " = E(- 1 ) n rtliK) 2n <0|z 2 "|0> (v-v') 2n 

n n - n \ln). 

(12) 

and we obtain the relations between the derivatives of the form factor and the 
moments : 

/W(0) = (-l) n p(m,) 2n < 0|^ 2 "|0 > (13) 
5 



or, by spherical symmetry, 

/ (n) (0) = ( _ 1) n_^L_ K)2 n /i2n (M) 

with 

M2n =< 0|r 2n |0 > . (15) 

2.2 Relation between the universal NR form factor and the 
NR limit of the Isgur-Wise function. 

Thus, the NR form factor F(P', P) is a function of the variable (v — v') 2 , while 
the relativistic Isgur-Wise function £(w) depends on w, where 

P P' 

w = v ■ v with v = — , v — —— (16) 

and (M,P), (M',P') are the masses and four- momenta of the initial and final 
mesons. 

There are several facts that ask for care in the identification between the NR 
form factor /((v — v') 2 ) and the Isgur-Wise function £(w) in its non-relativistic limit 
£nr(w). 

Let us leave aside for the moment the fact that the velocities (JTfij) differ from 
their non-relativistic limits (JHJ) by the binding energy that is neglected in the latter. 

Of course, /((v — v') 2 ) and £nr(w) cannot be generally identified, since w is not 
a function of (v — v') 2 : 



w = v ■ v' — v 1 + v /2 Vl + v 2 — v ■ v' 



1 

1+ 2 



v — V 



T- vT+v^-Vi + v 2 



(17) 



To relate w and (v — v') 2 one needs to choose a frame. The natural frame is the 
rest frame of the initial particle, i. e. v = 0. One has, in this frame : 



w 



1 + (v - v') 2 



1/2 



The relation between w and (v — v') 2 being non-linear, the relations between the 
derivatives relatively to w and to (v — v') 2 are complicated. As we show in Appendix 
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B, the derivative of order n relatively to w depends on the derivatives of order m < n 
relatively to (v — v') 2 and conversely. 

Another frame, the equal velocity frame (EVF) where the velocities are equal 
and opposite v' = — v, gives from (fTTj) : 

w = l+ l -{w-V) 2 . (19) 

Thus, in this frame, the relation between w and (v — v') 2 is linear, and the n-th 
derivative relatively to w is proportional to the n-th derivative relatively to (v — v') 2 . 
In the EVF, we obtain 

n n! 

= £nr(w) = (ip |exp [im q (y - v') ■ r q ) \ ip) 

OO 1 

= E(- 1 ) n 79^T K) 2 " <ok 2 "|o> (v-vf" 

n =o \ zn )- 

oo i 

= E(-!)" 2" T^y K) 2 " < 0|^|0 > (w - IT • (20) 

Therefore, in this frame one gets the relation between the derivatives of £,nr{w) and 
the moments : 

= (-1)" 2" J^Y), K) 2 " < 0|r 2 "|0 > . (21) 

This relation is exact in this frame, and coincides with the leading term in the NR 
limit in all collinear frames, as the rest frame (see appendix B). 

Therefore, in a NR expansion for the light quark, one can only claim to obtain 
frame-independent results for the derivatives of ^nr{w) in the leading NR order for 
the moments. From now on we will then rely only on the relation (|21|). 

However, v, v' in relation (|20|) are not identical to their non-relativistic limits 
(jOJ). This fact does not invalidate the relations given above, since w in the EVF 
(I19|) differs from its NR expression in terms of NR velocities by terms of order e/rriQ 
where e is the binding energy. This can be summarized by expanding w in terms of 
the NR velocities Q and binding energies : 

V 

2 \m q + m Q m q + m Q > 



1 / P P' 

w = 1 + I + Sub leading terms in velocities + O (e/ijiq 



(22) 



Therefore, since these differences are sub leading we can make the identification 

UM=/((v-v'f) (23) 



where wnr is the leading term of w given by the preceding expansion, and relation 
(|21|) holds indeed for the leading terms. 

From eq. (|21|) one may be surprised that the NR expansion leads to a result for 
£( n )(l) increasing with m q R ~ (t>/c) -1 , seemingly at odds with the notion of a NR 
expansion. In fact, the two last formulas in expression (|2T)|) show that the form factor 
is expanded in a series of powers (m q ) 2n < 0|z 2n |0 > [(v — v')/c] 2n ~ (m q R) 2n [(\ — 
V)/c] 2n where R is the bound state radius R ~ (j>int)~ l , Vint being the internal 
quark momentum. On the other hand v — v' is the external center-of-mass velocity 
transfer. Therefore, since m q R ~ m q /p int = [(t>/c) int ] _1 , the form factor is expanded 
in powers of the type [(v /c) ini ] _2n [(t> /c) ext ] 2n . The negative order in [(v /c) ini ] _2n is 
therefore compensated by the corresponding positive powers of [(v/c) ex t\ 2n . The 
derivatives £^(1) in equation (f2*T]) are given by the inverse powers [(v/c)i nt ]~ 2n , 
that multiply in (j2*Uj) [(v — v')/c] 2n = 2 n (w — l) n . In a general frame, the coefficient 
of (v — v') 2ri or (w — l) n will contain an expansion in powers of [(f /c)^ n^ ] _2n+m 
(m > 0). However, the subleading terms (m > 0) are frame-dependent, as shown in 
Appendix B. 

3 Constraints on the moments. 

Let us define the moments 

H n = < 0|r n |0 > (24) 
and consider the even moments fi2n, related to < 0|z 2n |0 > from rotational invariance 

ti 2n = (2n + 1) < 0|z 2n |0 > . (25) 

We will now formulate necessary constraints on the /i2n resulting from the fact that 
they are indeed moments, i.e. that there exists a function <p{r) such that 



poo 

H%n= / r 2n W)\ 2 dr 
Jo 
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(26) 



It turns out that these conditions are sufficient, but this is only proved in Appendix 
C, implying that the constraints are optimal. 

(a) A necessary condition is that for any non-zero polynomial P 



n roc n 

P(r 2 ) = W 2i > (r > 0) / P(r 2 ) \<p(r) \ 2 dr = £ a if x 2i > . (27) 

i=0 i=0 

From this condition, that is not very explicit, one deduces the following condi- 
tions (b), (c) and (d), that are equivalent. Condition (d) is explicit. 



(b) For any n > and non- vanishing do, • • • a n one has 



f*2i+2j > and ( a i)* a j ^2i+2j+2 > . (2£ 

i,j=0 i,j=0 



One demonstrates (b) from (a) by considering the polynomials P(r 2 ) = | a i r 



2i\2 



i=0 



lb 

and P(r 2 ) = r 2 \ J2 ciir 2% \ 2 . Conversely, (a) results from (b). 



i=0 



(c) For any n > 0, the matrices (fi2i+2j)o<i,j<n and {^2i+2j+2)o<i,j<n are positive 
definite. 

This condition is just a rephrasing of condition (b). 



(d) For any n > 0, one has 



det 
det 



(/*M+2;)o< i>3 -<J > (29) 

(/i2i+2i+2) < iii < n ] > . (30) 



To obtain (d) from (c) it is enough to note that a positive definite matrix has 
strictly positive eigenvalues, and that the determinant is the product of its eigen- 
values. 

Let us first write the determinants (J29j) and (J3U|) for the lower values of n, namely 

/i 2 > (31) 



det ( 1 ^ \ > (32) 

\^2 /i4 / 
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det 


V/^4 




> 






( 1 


/'2 






det 




//,! 


^6 


> 




\/i 4 


/'6 


/i 8 / 








/i 4 




\ 


det 


fJ>4, 


/'6 


/'s 


> 




V //« 




//in y 


/ 




f 1 


/'2 


/l.l 


A*6 \ 


det 


^2 


/i 4 


flf, 






/i 4 


/'6 


/'s 








//•8 


^10 


Vl2/ 



> 



(33) 



(34) 



(35) 



(36) 



where ()32|1. (J53J) and belong to the class of positivity conditions (|29|h and ()31|). 
and dHHJ) to the class flUJ). 
From (j^H) and (JSSj) we find, respectively 



/i 6 > 



/'2 



(37) 
(38) 



To get the constraint on /zg from (|34|) in terms of positive definite quantities, we 
make use of the following identities, 



//4 det 



/ 1 /X 2 A*4\ 

A*2 A*4 ^6 
V /i4 ^6 At8 / 



det 1 * * " "» - det 



.^2 A*4, 



. A*6 Ms . 



-1 2 



. A*4 A*6 / . 



(39) 



det 



/ 1 /i 2 /i 4 \ 

A*2 ^4 ^6 
V A*4 A*6 ^8 / 
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det 



/ A*4 1^2 \ 

(let I J det 



.M2 



1 /i 4 
/I4 /is 



det 



M2 



(40) 



/i 2 1 yU 4 
V/i 6 yU 4 fl 8 J 

that follows from the general identity among determinants of the Appendix A. We 
find : 



/ij - 2/x 2 /i 4 /x 6 + n\ _ (/i 2 /i 6 - f4 

> o — 7 — H 



(/i6 " ^ + (41) 



^4 — 1^2 ^4 (^4 — ^|) A*4 A*4 — A*2 

where the first equality follows from (|34j) and (|39|) and the second from (|34j) and 
(jUJ). 

To proceed in the same way with the 10-th moment, we make use of the inequality 
(135)) and the relations among determinants (A.8)-(A.10). 
The inequality ()35|). together with (A.8)-(A.10) yields 



A*io > 



A*2A*6 - A*4 



+ — 



A*6 (/^2yW6 - Ml) y"6 ^2 (^2/^6 - A*i) A*2 



+ (42) 



Things become more complicated for higher moments, but the method proceeds in 
the same way. 



4 Bounds on the derivatives. 

Let us summarize the inequalities among the moments deduced in the previous 

section. We adopt in (fTTj) and (jH)) the expressions given by the last equalities in 
the r.h.s. This will be instructive, as it will become clear below. We have obtained 

/i 2 > (43) 

M4 > (44) 

Me > - (45) 

M2 

(/i 6 - /! 2 /i4) 2 / 

yu 8 > 2 ^ A*4 (46) 

/i 4 — /X 2 

(/i 2 /X 8 ~ /i4/i 6 ) 2 , Ml /.-v 

/iio > : ay + — • (47) 

Ai2 (/U2/i6 - Ml) ^2 
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From the frame-independent relation between moments and derivatives obtained 
from (gTJ-flSS) ( cf - Appendices B and D) : 



[-l) n T 



\ 2n 



^2r 



(48) 



and from 
derivatives 



' NR ^' y ~' ~ (2n+l)! 
■ (|47|) . we obtain respectively the following inequalities among the 



(49) 



(50) 



-d 3 kD > ~ 



d 2 ki) 



■>NR\ 



(51) 



d 4 ki) > 



d 3 ki) + f a&a) d z ki) 



(i) 



.(2) 



(2) / 
AW 



5 r d 2 ki)l 2 



(5) 
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IT 



(1) 



t(2) 



1 N 5 SjVj; 

NRK L ) 9 



3£m 



t( 3 ) i 5 

— SJ 1 



£ (2) (1) 



21 

35 
99 



(3) / 
7VM 



(1) / 
NR\ 



(52) 



(53) 



siVi^l) + 7 

Importantly, we observe that the l.h.s. and the r.h.s. of all the inequalities 
(J53|) scale in the same way in the parameter R 2 m 2 , where m q is the light quark mass 
and R the bound state radius, since the derivatives f rom (EH), scale like 

(R 2 m 2 ) n . The inequalities (|4*9*|) and (JHUj) are the non-relativistic limit of the bounds 
of the true IW function © and (J3J). 

Comparing with the results of the method used in Appendix D, we have found 
stronger results for all the derivatives. 

For the even derivatives, we find for 2n > 2 a new term that strengthens the 
lower bound. This can be seen from the bounds flU and (D.9) for £$ R (l). We 
have found a new term besides the term proportional to [Cj^k-Q] 2 - However, for the 

(2) 

curvature ^(1), relevant for the non-relativistic limit of the curvature of the true 
IW function £(u>), we find the same bound (@J) as with the former simpler method 
of Appendix D. 
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As for the odd derivatives, the trivial bound (D.10) has been changed in a very 
substantial way, since we find the lower bounds (j51j) and (j53jl . 

Finally, let us emphasize that the lower bounds (}4*9"|) - (j5lH) are optimal. The 
optimality of these bounds is demonstrated in the Appendix C. 



5 Some illustrations. 

For the sake of a simple illustration of the bounds, let us consider the harmonic 
oscillator in the equal velocity frame : 



e N ° R H = exp [-(w - l)m 2 q R 2 \ (54) 

where the bound state radius R is normalized in a convenient way to have this simple 
expression. The n-th derivative reads 

dt(l) = (-l) n H^ 2 ) n • (55) 

Then, the bound (J49|) reads simply m 2 R 2 > and (J5U)) - (j53j) will become for the n-th 
derivative 

n = 2 1 > - 
5 

n = 3 1 > \ 

7 

1 55 
n = 4 1 > — 

63 
91 

n = 5 1 > — . (56) 

99 v ; 

Interestingly, we find that the bounds become better and better as we consider 
higher derivatives. For the 5-th derivative the bound is already very strict. 

However, it is not granted that these features will remain for more realistic 
potentials. Therefore, it can be useful to examine another simple potential, although 
not confining, namely the Coulomb potential. In this case we have a dipole form 
factor 



f N Oulo mb{w) = (57) 

l + {w- l)mlR 2 
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The derivatives read 



= (» + !)! (-1)" HR 2 T (58) 
and the inequalities (|5U |) -(|53 )l give, respectively 









2 


n 


= 2 


1 > 










5 






1 > 


15 


n 


= 3 


28 






1 > 


263 


n 


= 4 










378 






1 > 


3626 


n 


= 5 








4719 



(59) 



These inequalities are somewhat less strict than in the harmonic oscillator case 
but here also they improve for higher derivatives. 

We can expect that in the case of a realistic phenomenological Qq potential, with 
a confining and a short distance parts, the situation will be in between the harmonic 
oscillator and the Coulomb potentials. 



6 Generalization of the method to QCD. 

We have obtained lower bounds on the derivatives at zero recoil of the non- 
relativistic Isgur-Wise function £nr(w), i.e. the IW function with a NR light quark. 
Our main motivation has been to find the leading term in a NR expansion of the 
derivatives at zero recoil of the true IW function £(w) that should be obtained in 
the heavy quark limit of QCD. The parameter in this expansion is (v 2 /c 2 )i n t or, 
equivalently 1/R 2 m 2 , where m q is the light quark mass and R is the bound state 
radius. In previous work [3] we did obtain in the heavy quark limit of QCD such an 
expansion for the slope and the curvature, inequalities (J2J) and 

3 



-e (1) (l)>, 
-( 2 )m >s 1 



^(1)> -|-4^(l) + 3^ ij (l)J | . (60) 

Since — £^(1) and ^ 2 \l), scale respectively like R 2 m 2 and (R 2 m 2 ) 2 , in the NR 
limit these inequalities become respectively (|4T?|) and (f57Ij) . The inequalities (f6T)j) 
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contain terms, specific to QCD in the heavy quark limit, that are subleading in a 
NR expansion. 

Our aim would be, in the long run, to obtain bounds for the n-th derivative of 
the IW function in the heavy quark limit of QCD that must contain the subleading 
terms in a NR expansion. We know that in the strict NR limit we must recover the 
bounds (jHH)-(J5Sl) obtained in the present paper. 

To obtain these bounds in QCD we could try the method of [3\ in a systematic 
way, that uses sum rules for the non-forward amplitude, relating a sum over inter- 
mediate states and the OPE, that depends on three variables Wi — Vi -v ', Wf — Vj ■ v', 
Wif = Vi ■ Vf, that lie in a certain domain jlj. Differentiating the SR relatively to 
(wi,Wf,Wif) and going to the frontier of the domain one gets relations that allow 
to obtain ()60|) [3 . This method can be pursued further and obtain bounds for the 
higher derivatives [Zj. However, the obtained bounds, in their NR limit, are weaker 
than dSU-dSnj)- 

We have developed here a more powerful method, based on the positivity of 
matrices of moments of the ground state wave function, that allows to go further for 
the derivatives n > 2 in the NR limit. To generalize the present method to QCD 
in the heavy quark limit one should investigate whether the derivatives £^(1) can 
be expressed in terms of positive definite quantities that are true moments as in 
the non-relativistic expression (|26|) . Then, one could draw the consequences that 
follow from the positivity of the relevant matrices. A step in this direction is the 
conjecture that, at least in the meson case £Q, all SR in the heavy quark limit of 
QCD are satisfied in the Bakamjian-Thomas (BT) class of relativistic quark models 
[Oj. We have realized this in practice for the lower derivatives £( n )(l), n = 1,2,3. 
These models are relativistic for the states and also for the current matrix elements 
in the heavy quark limit, exhibiting Isgur-Wise scaling. One can hope to start from 
the NR quark model and go to BT models, and from those to the heavy quark limit 
QCD. 

Another, more direct way to proceed to the heavy quark limit of QCD is to 
start from the sum rules obtained in pQ-jH], realizing that one can obtain the NR 
bounds of the present paper from the equivalent sum rules of the non-relativistic 
limit. Indeed, in the NR limit we have a SR of the form 
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E fn,n»(k)fn»,n>(U) = /n,»'(k + k') . (61) 
n" 

that follows, very simply, from 

f n ^{k) = <n\e i ^\n'> . (62) 

In QCD in the heavy quark limit we have sum rules of the same form ([61)1 . but 
without the explicit expression ()62|) . However, to derive the inequalities of this 
paper, (JoTj) is sufficient. Indeed, from (JoTj) one gets 

5>c*/ ,o (ki - k,) = EE c ^/o,n(k,)/„, (-k J ) = E | E c ^o,n(k,)| 2 > . (63) 
From this relation one can infer, for any function y?(k), 

J dk dk' <f(k')* / ,o(k - k') <^(k) > (64) 
and therefore, for the Fourier transform 

|rfr|^(r)| 2 / , (r)>0 (65) 

and hence 

/o,o(r) > . (66) 

These are the conditions that we need to obtain constraints on the moments and 
hence bounds on the derivatives of the form factor /o,o(k — k'), because writing the 
form factor in terms of its Fourier transform 

/ , (k-k') = Jdr / 00 (r) e^"^ (67) 

and taking into account that /o,o(k — k') must be an even function, we obtain (Oz 
is defined along the momentum transfer k — k') : 

/oo(k-k') = E(-l) n t-^t (k-k') 2n (68) 
J0 '° K ; V 2n ! v ; 2n + l y J 



16 



i.e. an expansion of the form factor in terms of moments of the form (|20|). with the 
identification 



/i 2 n = < 0|r 2 "|0 > = Jdx / ,o(r) r 2n (69) 
since the positivity condition /o,o( r ) > (JBTIj) holds and, from (jBTf) : 

J /o,„(r)dr = 1 . (70) 

We do recover essentially the previous results ()43j) - (j47j) usm<? onZy t/ie sum ru/es 
f!61|) . If /o,o( r ) is a function, it can be seen as the square of a wave function, and all 
the strict inequalities of the NR type (J4"!^ - (J4"7]) would follow. However, this is not 
implied by the SR, and weaker results could follow, namely the inequalities may not 
be strict. For example, one could have a distribution like 

/o,o(r) = ^(|r|-r ) (71) 
which is not the square of a wave function, that would imply 

Urn = (72) 

and, for example, the strict inequality (J44j) would become an equality. For example, 
in the true QCD case, the lower bound (J2J) for p 2 could become an equality. By the 
way, this would correspond to the so-called BPS approximation [T2"] . 

Our strategy will then be to start from the SR in the heavy quark limit of QCD 
that are equivalent to the NR ones (|6*Tj) . and proceed along the same lines. We can 
presume that the method will give the optimal bounds for the derivatives of the true 
Isgur-Wise function £(u>). 

7 Conclusion. 

To conclude, we have obtained the best possible general bounds on the deriva- 
tives of the Isgur-Wise function (.nr^), i.e. considering the light quark as non- 
relativistic, in terms of lower derivatives. These bounds must be the non-relativistic 
limit of the bounds on the derivatives of the true Isgur-Wise function C,(w), and 
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constitute a guideline in the derivation of the latter. Moreover, we argue that the 
method developed here, that exploits the positivity of matrices of moments can be 
generalized, starting from SR in the heavy quark limit of QCD, to obtain the best 
bounds on all the derivatives of £(w). 
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Appendix A. An identity between determinants. 



In Section 3 we have made use of the identity among determinants 

det [(o^) l<i,j<n] det [{o>ij)2<i,j<n-l\ — det [{a>ij)l<i,j<n-l\ det \{flij)2<i,j<n\ 
— det [(aij)l<i< n -l,2<j<n] det [((lij)2<i<n,l<j<n-l] 

or, in a more readable way : 



(A.l) 



det 



/ ai,i ai j2 

&2,1 ^2,2 

V a n,l &n,2 
/ «2,2 



01, n -l <3l,ra \ 

02, n-l a 2,n 



= det 



det 



a n-l,2 
V On,2 

a n-2,2 ' 
V a n-1,2 ' 



0"n—l,n—l @"n—l,r 



/ a-2,1 



det 



0-2,n-l \ 
Cln-l,n-l ) 



Va n -1,2 • 

«2,n-l «2,n \ / " " " «l,n-2 dl,n-l \ 



Q>n— l,n— 1 Q"n—l,n 
&n,n—l dn,n 



det 



a n-2,l 
V On- 1,1 



Ql,n-1 &l,n \ / 02,1 



a n-2,n-l a n-2,n 
a n—\,n—\ ®"n— l,n / 



det 



Qn-1,1 
V On,l 



a n-2,2 a n-2,n-l 
O-n-1,2 &n-\,n-\ } 
0>2,n-2 Q2,n-1 \ 



a n-l,n-2 a n-l,n-l 
^n,n— 2 Q>n,n— 1 / 



(A.2) 



To demonstrate this relation, let us introduce the column vectors 



i=i 



(A.3) 



where the Ojj are the elements of the matrices (A.l) or (A.2). Multipying (A.2) by 
(ei A • • • A e n ) <g> (ei A • • • A e„), this formula writes 

(xi A x 2 A • • • A a; n _i A x n ) (8) (ei A x 2 A • • • A x n -i A e n ) 
= (xi A x 2 A • • • A x n -i A e„) <g> (ei A x 2 A • • • A x n _i A x n ) 
- (e n A x 2 A • • • A x n _i Ai n )®(iiAi 2 A---A x n _i A ei) . (A.4) 
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Assuming that the vectors independent, one can expand e\ and e 



ei = aiXi + • • • + a n x n 
The l.h.s. of (A. 4) becomes 



e n = faxx H h /3„a; r . 



(A.5) 



{x\ A x 2 A • • • A x n _i A x n ) ® (ei A rc 2 A • • • A x n _i A e n ) 
= (xi A x 2 A • • • A x n _i A x„) <g> ((aiXi + a n i n ) Ax 2 A---A x n _i A (ftzi + /3 n x n )) 
= - (xi A x 2 A • • • A x n _i A x„) <g> (xi A x 2 A • • • A x n -i A x n ) (A.6) 

while the terms in the r.h.s. become 

(x 1 A x 2 A • • • A x n _i A e n ) ® (ei A i 2 A ■ ■ ■ A x n _i A x n ) 

= ai/3 n (xi A x 2 A • • • A x n _i A x„) ® (ii A i 2 A ■ • ■ A x n _i A x n ) 



(e„ A x 2 A • • • A x n -i /\ x n ) ® (xi /\ x 2 /\ ■ ■ ■ /\ x n -i A ei) 
= (xi A x 2 A • • • A x n _i Ax n )®(iiAx 2 A---A rc n _i A x n ) . (A.7) 

The identity is therefore demonstrated if the vectors X\ i ' ' ' vCfi £1X6 independent. 

If the vectors are dependent one can show that both the l.h.s. and the r.h.s. 
vanish identically. 

In particular, we make use in Section 3 of the following identities : 



det " 4 " 6 de* 

V A*6 A*8 / 



det 



//i 4 fiQ // 8 \ 

V/i 8 /iio A*12 / 



Z' 1 pL 2 H4 fX 6 \ 

A*2 A*4 A*6 ^8 

A*4 ^6 ^8 A'lO 

\/i 6 /i 8 /iio /i i2 / 

/ 1 /i 2 /i 4 \ 



det 



A*2 A«4 A*6 



\ 


/ 


/ 


det 

V 



A*4 A*6 A*8 



(I /i 4 \ 
det det 

\A*4 A*8 / 



/ /i4 A*2 A*6 /i8 \ 

/X 2 1 /i 4 /i 6 

A*6 A*4 ^8 A*10 

\fi 8 /i 6 /iio A*i2 y 
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(A.8) 



/ A*4 A*2 A*6 \ 

det /i 2 1 /i4 
V yU6 A*4 A*8 / 

det 



det 



/ 1 // 4 
A*4 A*8 A*io 



\ 




/ 


det 

. V 



fi 2 1 Pa\ 

l-l>6 A*4 /^8 



/i4 /i 2 
.A*2 1 . 



det 



/ ^8 ^6 A*4 A*10 \ 
A*6 ^4 ^2 A*8 

V/ilO yU8 A*6 Ail2 / 



det 



/ /^8 ^6 A*4 \ 

/i6 /i 4 /i2 det 

V/i 4 /i 2 i / 



A*2 1 A*6 



\ 


/ 


/ 


det 



A*4 A*12 



(A.9) 



(A.10) 



Appendix B. Frame-dependence of the subleading 
moments. 

Let us begin from the general expression (jlOj) . (fTHj) in the initial hadron rest 
frame, v = 0, without neglecting powers < 0|2; 2p |0 > (w — l) p (p < m). From the 
relation 

(v - v') 2 = w 2 - 1 = jr 2 m_< ( m ) {w - l) m+l (B.l) 

i=0 V 1 ) 

one obtains the following relation between derivatives and moments : 

4t(i) = ± (-ir 2 2 — w!t "; T K) 2m < o|, 2 -|o > . 

m t^/2 ( 2m ) ! ( n - m ) ! ( 2m ~ n > ■ 

(B.2) 

conversely, to obtain the moments in terms of the derivatives, it is enough to replace 
the variable w by 

y = w 2 - 1 . (B.3) 

Then 

&v*H = £W(l+2/) 1/2 ] (B.4) 
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can be expanded in powers of y = w 2 — 1 : 



W(i + y) 1/2 ] = E 3 E cn, m y 



n=0 



n! 



OO 1 

= Y(-i) m —— 

h [ } (2m)! 



(m g ) 2m < 0|^ m |0 > y m . 



2m I 



(B.5) 



From this expression we can read the expression of the moments in terms of the 
derivatives 



m i 

< 0|^ m |0 >= (-l) m (2m)! K)- 2m E ~, c ^ dt(l) 
The coefficients c„ irrt are defined by 



n=0 



i+y) 1/2 -i] n = Y,cn, m y r 



and therefore given by 



n 

"n,m = E( _ 1) 
i=0 



'raN /i/2 N 



i I \ m 



This sum can be calculated and gives : 



(-1) 



m—n r\—2m+n 



2m — n — 1\ (2m — n — 1 
, m — n J \ m — 77 — 1 , 



Explicitly, one has 



(B.6) 



(B.7) 



(B.8) 



(B.9) 



Co,m — ^m,0 
( _ 1)m -n 2 -2m+n n (2777 - 77 - 1) 

m\(m — 77)! 



(B.10) 



where the second relation holds except for 77 = m = 0, since c 0j o = 1- 

Gathering relations (B.6)-(B.10) we obtain the final relation giving the moments 
in terms of the derivatives 



< |^ 2m |0 >= KP m 5 mfi 6v R (l) 



2m 



. f_1 \n 0-2m+n (2m) ! (2m - 77 - 1) ! {n) 

+ ti m\{m - n)\{n - l)\ W j " 

The relations (B.2) and (B.ll) are the main results of this section. 



(B.11) 
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Explicitly, one obtains, for the lower derivatives and moments : 



^ fl (l) = <0|l|0> 



1 



-m\ < 0\z 2 \0 > 



= 3 K < 0\z 4 \0 > - m\ < 0\z 2 \0 > 

d'kl) = ~ rn* < 0|/|0 > + m\ < 0|/|0 > 

1 2 

ZnrO) = 77^ m a < °k 8 |° > -7 m a < I^ 6 |° > + m t < O^l > 
105 9 5 9 9 



d^l) = ~ rnf < 0|/°|0 > + A m 8 , < 0|/|0 > -m\ < 0|/|0 > (B.12) 
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<0|l|0>=6v«(l) 
< 0|z 2 |0 >= -rn' 2 
<0|/|0>=3m- 4 - &*(!)] 

-d 3 ki) + 3 - 3 CWl) 

< 0|/|0 >= 105 m- 8 [dkl) " 6 d 3 i(l) + 15 - 15 



< Q|/|0 >= 15 nig 6 



< 0\z 10 \0 >= 945 m. 



-10 



-d 5 ki)+io6 d 4 ki)-45 dki) 

+ 105 - 105 & fl (l)] . (B.13) 

We observe that in expressions (B.2) and (B.ll) there is a leading term in the 
non-relativistic expansion and in (B.12) and (B.13), the first term in the expansion 
is the leading term. 

The purpose of the present detailed calculation is to point out that actually the 
subleading terms are frame- dependent and that therefore, one can only get infor- 
mation on the leading term in the non-relativistic expansion, the unique term that 
appears in the equal-velocity-frame used in Section 2. 

Moreover, there is continuity between the rest frame and the equal-velocity 
frame, since w can be expressed in terms of (v — v') 2 in all collinear frames 



av + /V = , a + fi = l 



(B.14) 



one gets indeed 
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w = yj\ + (1 - «)2( v - v') 2 yjl + a 2 (v - v') 2 + a(l - a)(v - v') 2 (B.15) 

with a = 1 in the rest frame and a = \ in the equal-velocity-frame. The first order 
in (v — v') 2 yields, independently of a, 

w 2* 1 + -(v- v') 2 (B.16) 
giving the leading order relation between derivatives of £nr(w) and moments (f2*Tf . 



Appendix C. Optimality of the constraints 

We have seen that the non-relativistic Isgur-Wise function 

1 



^H = <^ |e-^ (v '- v) - r |^o> 
has its derivatives at w = 1 related to the moments 



W = ^(v'-v) 2 + l 



(c.i; 



/i 2 n = < Mr 2n \^ > (C.2) 

by 

= (-1)" ^^yy mf ^ (C.3) 
and that these moments satisfy the following tower of constraints : 

det (jkK+2j) < ijj < n > (C.4) 

det (/i 2 i+2i+2)o< i , j < n > (C.5) 

which consists, for each n > 0, of a lower bound on // 2n depending on the moments 
/j,2k for < k < n — 1. 

In this appendix, we show that this cannot in general (for arbitrary wave func- 
tion) be improved. Namely that, given \x-ik for < k < n — 1 satisfying (C.4) and 
(C.5), the moment /i 2n can have any value larger than this lower bound. 
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To that goal, we forget (C.2) by now, and consider arbitrary numbers /i 2n (n = 
0,1,2, •••) satisfying the constraints (C.4) and (C.5). In this appendix we prove 
that, for any N > 0, there exists a wave function such that one has 

^ 2n = < |r 2n |4 W) > for all > k > N . (C.6) 

We shall not be able here to know if there is a wave function ip satisfying (C.2) for 
all n, but our more limited result (C.6) is enough to prove the point. 

To simplify notations, we introduce 

x = r 2 (C.7) 

as a variable taking positive values. 

1 - Introduce, in the vector space of polynomials in x, the linear form defined 
by the values fi2k on the monomials x k (which constitutes an algebraic basis of this 
vector space). It is given by 

n n 

< P > = a k A*2fc for P(x) = J2 a k xk - ( C - 8 ) 

k=0 k=0 

As a preliminary and crucial step, we have : 

2 - This linear form < P > is strictly positive. Namely, one has 

P(x) > for all x > 0, and P ^ < P > > . (C.9) 

To prove this, observe first that (C.4) and (C.5) imply that (^2i+2j)o<i,j<n arid 
(fi>2i+2j+2)o<i,j<n are positive definite matrices, or explicitely that one has 

n n 

a% a* fi 2 i+2j > , Oi a* fi 2 i+2j+2 > (CIO) 

i,j=0 i,j=0 

for any coefficients ao, • • • , a n not all vanishing, and using the definition (C.8) of 
< P >, these properties (CIO) translate into : 

<|Q| 2 >>0, < rr|Q| 2 > > (C.ll) 
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for any no n- vanishing polynomial Q. 

Then, any P satisfying P(x) > for all x > is a linear combination with 
positive coefficients of polynomials of the form \Q\ 2 or a;|<5| 2 . Indeed, considering 
the roots of P, we have 

P = Y[.{x + q) H.(x - c;.) 2 II fe \x - z k \ 2 (C.12) 

with Q > 0, c'j > 0, Im z k ^ 0, since complex roots z k occur in conjugate pairs, 
strictly positive roots c'j occur in event multiplicity (else a change of sign at x — cj-), 
and negative roots — q are arbitrary. 



3 - Next introduce a scalar product in the vector space of polynomials by 

< p\p> > = < p*p' > . (C.13) 

The scalar product properties are easily verified. Notably, the important fact that 
< P\P >= implies P = results from (C.4). 

We may then consider the orthogonal polynomials Po,Pi,P2, • ■ ■ with respect to 
this scalar product. The theory of orthogonal polynomials is classsical matter . 
They are usually considered with respect to a scalar product defined by a weighted 
integral, but their properties extends easily to the more general case needed here, 
where we do not know a priori if the scalar product (C.13) can be given by an integral. 



The polynomial p n has degree n, and we have : 



< Pn\Pn> >= d 2 n 5 n , n i . (C.14) 

It will be convenient for us to fix p n by taking the coefficient of x n to be 1. These 
polynomials can be computed recursively by the orthogonalisation Schmidt process : 

Pn = X n -Y j <M f > Vk (C.15) 

t^Q < Pk\Pk > 

where the automatic fact that < pu\Pk >7^ is essential. Also, since any polynomial 
of degree < n is a linear combination of po, ■ ■ ■ ,p n , we have the property : 
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< P\p n > = for any P of degree < n (C.16) 

which is of constant use in the following. Taking P = l,x, • ■ ■ , a;™ -1 , (C.16) gives 
a system of n linear equations for the n + 1 coefficients of p n , which, according to 
(C.4), can be solved uniquely up to a constant, and then Cramer's formulae give an 
explicit expression for p n . 

We are actually interested by the zeros of p n . 

4 - All the roots of p n are simple and strictly positive. 

In fact, let be the strictly positive roots of p n of odd multiplicity. We 

have only to show that m = n. If m < n, according to (C.16), we have 

< (x - Xi) ■ ■ ■ (X - X m )p n > = < (x - Xi) ■ ■ ■ (x - X m )\p n > = . (C17) 

However, the polynomial ( )p n has a constant sign for x > 0, and 

does not vanish identically. Therefore, according to (C.9), (C.17) and hence m < n 
cannot be. 

5 - We may now write explicit formulae for /i 2 k with < k < 2n — 1 : 

n 

H2k = Y, X i x i (0<k<2n- 1) (C.18) 
i=i 

where the roots of p ni and the coefficients A, are given by : 

i <E M > (C19) 

p' n {Xi) X-Xi 

To prove (C.18), notice that it amounts to : 

n 

< P> =J2 X i p ( x i) ( an y P of degree < 2n - 1) . (C.20) 

i=i 

Performing the Euclidean division of P by p n , we have 
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P = Q Pn + R (C.21) 

with degree Q < n and degree R < n. We may verify (C.20) separately for Qp n and 
for R. 

For P = Qp n , the left-hand side of (C.20) vanishes by (C.16), and the right-hand 
side vanishes because the the roots of p n . 

For P = R, we use the identity 

fl M = £^f=T ,iW (C - 22) 

which stems from the fact that both sides are polynomials of degree < n, that are 
equal at n points x = Xj. Then (C.20) is satisfied due to the choice (C.19) of the 
coefficients Aj. 

6 - Define 

n 

Pn (x) = £A; S(x- Xi ) . (C.23) 
i=i 

Then (C.18) writes 



poo 

V2k = / p n (x) x k dx (0 < fc < 2n - 1) (C.24) 
Jo 

to be compared with formula (C.6) to be proved, which writes 



roo 

P-2k 



roo 

/ PN {x) x k dx (0 < k < N) (C.25) 
Jo 

with 



p N (x) = 2ir ^ jtf^v^)! 2 • (C26) 

We have still the problem that p n (x) is not a function. 

The idea to solve this problem is to vary p± n -2-> keeping fixed p 2 k for < /c < 
2n — 2. The polynomial p n then depends on p^n-i as a parameter, and as well 
its zeros Xi{p± n -2) and the coefficients Xi(p4 n -2) defined by (C.19). Then formula 
(C.24) is lost for pi n -2, but remains valid for < k < 2n — 2, and in fact gives a 
whole family of formulae 
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V2k = / P«(/i4n-2, x) x k dx (0 < k < 2n - 2) (C.27) 
with a weight distribution 

n 

Pn(l^4n~2, x) = J2 X i (/^n-2) 6(x - X i (fi 4n - 2 )) (C.28) 
i=l 

depending on \i± n -2- We may then take the mean value of (C.27) over any interval 
[/4n-2> / i in-2] i n which the constraints are satisfied, obtaining 



^2k 



/ p n (x) x fe dx (0 < fc < In - 2) (C.29) 



with 



,(2) 



Pn{x) = -pj (J) — / (1) Pn(A*4n-2,aj)dA*4n-2 • (C.30) 

Now, p„(x) defined by (C.30) has a good chance to be a genuine function, because 
integrating a 5 distribution over a parameter usually gives a function. 

However, there is an obvious case in which this does not hold, namely when the 
point where the 8 distribution is concentrated does not depend on the parameter. So 
we still have to show that each zero of p n does vary with yU 4n _ 2 . Let us consider the 
orthogonal polynomials po,pi, • • • with respect to the new scalar product < P\P' > 
associated to new values fio,l<>2, • • • of the moments, with p, 2 k = P2k for < k < 
2n — 2, and /Lt„_2 ^ Pi n -2- Note that the new scalar product of two polynomials 
is the same as the original one when the sum of the degrees is < 2n — 2. It follows 
that pk — Pk for < k < n — 1, and also that 

< Pn\Pk >=< Pn\Pk >= for 0<k<n-2 . (C.31) 
Therefore, the expansion of p n over the pk writes : 

Pn=Pn + C p n -! . (C.32) 

And one has 0. Indeed, since < p n \p n -\ > = 0, one has 

C = ) (< Pn\Pn-l > ~< Pn\Pn-l >) 

< Pn-\\Pn-\ > V 1 
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so that 



(C.33) 



c = : (/i 4 „+ 2 - • (C.34) 

< Pn-l\Pn-l > 

The fact that a zero of p n cannot be a zero of p n now follows from (D.32) and 
the fact that a zero of p n cannot be a zero of p n -i- 

This last point is a well known property of orthogonal polynomials, which can 
be proved directly as follows. Assume that a zero x, of p n is also a zero of p n -\- 
Then we have 

< Pn-ll-Z 5 - > = < > = (C.35) 

by (C.16). On the other hand, writing = ax 71 ' 1 + •■■■, where a 7^ 0, we have 
again by (C.16) : 

<Pn-l\ —> = a < Pn^x" 1 ' 1 > = a<p n ^ l \p n ^i> (C.36) 

X — Xi 

which cannot vanish, contradicting (C.35). 

We are now in position to complete the proof of (C.6). Indeed, using the implicit 
functions theorem, one can infer from (C.32) and (C.34) that, for a small enough 
interval [/4n-2> /4n-2]> eacn function Xj(/i4„_2) is a diffeomorphism of this interval to 
an interval [xf\ xf\ in x. Then introducing the reciprocal function Xi — > /ij(xj) of 
HAn-2 a ; i(/ i 4n-2), the integral of a 5 function is computed by changing the variable 
of integration fi4n-2 to X; = Xi(/i 4n _ 2 ) : 



(2) (2) 
/ 4 " 2 Ai(/i 4 n-2) 5 (x - Xi {Hin-2)) dfl 4n „ 2 = ' X(Hi(Xi)) 6(x - X t ) |//-(Xj)]dXj 



= Xu(i) ^,(1) A< (/ii(x)) |/4(x)| (C.37) 

where xi is the characteristic function of an interval J, namely Xi{ x ) = 1 for x £ I 
and X/( ;r ) — for x £ I. Then (C.30) gives 



Pn{x) = — } (jj— X! X^w^/x) Ai(/ii(x) |/i-(x)| (C.38) 

pin-2 ~ / i 4n-2 i=l 



1 ' 1 
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and this is a genuine positive function, which can therefore be written as the square 
of a wave function. 



Appendix D. Explicit lower limits for the even deriva- 
tives. 

In this Appendix we generalize to all even derivatives ^ 2n '(l) the proof of the 
bound (|5U|) that we have given in ref. jSj- 
From expression (|21|) : 

41(1) = 2 m K) 2m < 0|^|0 > (D.l) 

using rotational invariance we obtain 



and from 



tFAO) = ("l) m T^y K) 9 — JT < 0|r 2m |0 > (D.2) 
(2m) ! 2m + 1 



< 0|r 2m |0 > = |< 0|r m |0 >| 2 + £ |< n|r m |0 >| 2 . (D.3) 



Using again spherical symmetry 



< 0|r 2m |0 > 



one obtains 



and therefore 



(m + l) 2 |<0|,2 m |0>| 2 + (m + l) 2 \<n\z m \0>\ 2 

nj^0,rad 



ml 



(2m)! 



m„ 



i2m 



(m + l) 2 
2m + 1 



< 0\z 2m \0 > 



+ \<n\z m \o>\ 

nj^0,rad 



(DA) 



(D.5) 



Wm. m! ^ ^2 m + i <0 i 2 m| 0> |2 



:-irc^(i)> 



(D.6) 



(2m)! v '" y/ 2m + 1 
This expression demonstrates that £nr(w) is an alternate series in powers of 
(w-1). 
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Assuming m to be even, m = 2n, one gets 



d 2 £(i) > 



(2n) 



(2n + 1) 



< 0k 2n |0 > 



(D.7) 



(4n)! v "* / 4n + l 
The moment < 0|z 2n |0 > can be expressed in terms of £^(1) though (D.l), giving 
finally 



[(2n)!] 3 (2n + l) 



[n!] 2 (4n)! An + 1 
We obtain, for the lower values of n, 



(n > 0) 



(D.8) 



n 



n 



r( 4 ) m >s 5 fc( 2 ) ' ll2 



1 > 



21 



(D.9) 



The formula (D.8) generalizes the result (|5U|) to all even derivatives. 

For the odd derivatives one gets, with the present method, from (D.6), the weaker 
result 

-dS +1) (l)>0 (n>0). (D.10) 

We see that we had obtained in Sections 3 and 4 a much stronger result for 
fjfi(l) than (D.9), and non trivial results for —£jv^(l) an d ~£7vj?(l)- However, we 
have obtained here an explicit lower bound for £^(1) (D.8). 
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